Apparatus and Method for Accurate Energy Device State-of-Health (SoH) Monitoring

ABSTRACT

A device and associated testing method for testing state of health of an energy device, comprising: applying electrical excitations to the energy device at a predetermined electrical excitation frequency ω e ; applying mechanical excitations to the energy device at a predetermined mechanical excitation frequency ω m ; measuring an electrically-induced phase difference Δθ e  (ω e ) between voltage (V) and current (I) within the energy device from applying the electrical excitations; measuring a mechanically-induced phase difference Δθ e  (ω m ) between voltage (V) and current (I) within the energy device from applying the mechanical excitations; and deducing the state-of-health of the energy device by comparing the electrically-induced phase difference Δθ e  (ω e ) with the mechanically-induced phase difference Δθ e  (ω m ).

BACKGROUND OF THE INVENTION

Electrochemical impedance spectroscopy (EIS) has been used in studies of electrode and plate behavior during charging and discharging for years. The impedance response of the battery, or more broadly, any electrochemical energy storage and/or conversion device (including certain types of fuel cell technologies, solar cells, and even certain types of capacitor technologies), depends on the measurement frequency and the “state” of the energy device at the time it is tested. For the specific case of a battery, this method has been reported to be affected, to varying degrees, by many fundamental battery parameters, including battery design, manufacturing tolerances, aging, temperature, and state-of-charge (SoC). The goal of EIS State-of-Health (SoH) monitoring techniques is to extract as much information as possible about the electrochemical health of the battery. Based on the variability introduced by the parameters noted above this monitoring technique has fallen short in the market place. The evolution of the technology followed a very predictable path and this application discloses the next generation of this technology which is a revolutionary break-through in the science of EIS SoH battery monitoring.

This disclosure will elucidate the novelty, inventive step and industrial viability of a method and device for more accurately measuring the SoH of an energy device by illustrating the drawbacks of the current SoH measurement and prediction techniques on lead acid battery chemistry, and showing how using ultrasonic energy as part of the measurement process overcomes these deficiencies. However, it should be understood that the disclosures here can be applied to any electrochemical energy storage and/or energy conversion devices (generally referred to simple as “energy devices”), including but not limited to, certain types of fuel cells, solar cells and capacitor technologies. In the broad understanding of the method and device disclosed here, it applies to any electrochemical energy device in which there are ionic species kinetic limitations and/or deficiencies. Certain fuel cells, solar cells and capacitors also have known ionic kinetic deficiencies which impact their operating characteristics and performance and it should be understood that this method and device disclosed encompasses these type of devices as well.

A brief background is necessary in order to elucidate the novelty of the invention being presented. In recent years there has been considerable activity and debate regarding the use of internal “impedance” characteristics as a battery condition measurement. The interest reflects the desire for simple electronic means to replace discharge testing as a practical determination of residual battery capacity, particularly given the increased usage of valve-regulated lead-acid (VRLA) batteries (for the lead acid chemistry example). It is widely accepted that the only way to truly know the true health of a battery is to periodically conduct test discharges. Once the test discharge is complete, and passes, this is a ‘good’ indicator that the battery is operating as designed and it provides peace-of-mind that the battery will perform its intended function when called upon. Unfortunately, conducting a test discharge is costly and time consuming. More importantly, if a test discharge is run on a battery pack too often it will ultimately reduce the battery's operating life, and the subsequent replacement will prove to be extremely costly. Batteries do not respond well to being fully discharged and it directly takes away from the battery life.

When conducting a test discharge it is often necessary to take the battery off-line, i.e. disconnect from its intended function, and therefore alternative measures for back-up power need to be in place. Typically, a back-up battery pack, or other means of power, such as diesel generator, fuel-cell, etc., is added to the system design to minimize the impact of this necessary operational restriction, however, this brute force attempt adds significant up front capital and maintenance cost to the project. The length of a test discharge for a primary battery back-up plant is dependent on the rate at which it is discharged. Typically, this time is determined by how much of a load is on the system and how long it may be needed in an emergency situation. As a point of reference, some applications require the battery pack to support critical loads for a period up to ˜72 Hours. Therefore, the minimum time the battery pack will be off-line is approximately 3 days, while the battery is being discharged. However, after the 3-day test discharge is complete the battery is still not available for its intended function. The battery must be charged back up and this is also a time consuming process. If the battery is charged too hard or too fast it will take away from the battery's operating life, and again the replacement costs are significant. Thus, a general rule is that the battery should be charged for at least the length of time it was discharged, however, battery manufacturers would prefer a much longer charge time (on the order of 30 days for very large lead-acid battery packs). Therefore, the minimum time the critical battery back-up source in this example will be off-line, excluding set-up and testing time, is approximately 6 days. However, this aggressive charging algorithm is on the edge of what battery manufacturers are willing to warrantee, so prudent engineering judgment would dictate using longer charge times, thus increasing the down time of a critical energy source.

A purpose of this novel, inventive and industrious mechanical excitation technique is to remove this inherent uncertainty in the electronic measurement technique, thereby obviating the need to conduct frequent and costly test discharges. Conducting periodic test discharges of the plants critical battery back-up not only places sever operational restrictions on the plant, but, even more importantly it reduces the operating life of the battery pack and results in premature battery pack replacement. If the utility tries to minimize the length of time that this necessary operational restriction is in place, by charging the battery at a faster than recommended rate, than this will ultimately result in unwanted and premature costly battery pack replacement. Thus, the peace of mind that comes with knowing that your critical battery back-up source is ready for its intended function comes at a great cost. The utility has to weigh the risk of minimizing operational restrictions against the looming cost of replacing its battery pack. As a result, many companies are attempting to determine the battery health without test discharging. The method disclosed herein will offer a piece of mind to the plant owners that there battery pack will perform its intended function without the need to perform costly, restrictive and timely test discharges.

SUMMARY OF THE INVENTION

A device and associated testing method for testing the State-of-Health (SoH) of an electrochemical energy storage and/or energy conversion device (“energy device”) and ascertaining addition electrochemical health information that would not be available without said method, comprising: applying electrical excitations to the energy device at a predetermined electrical excitation frequency ω_(e); applying mechanical excitations to the energy device at a predetermined mechanical excitation frequency ω_(m); measuring an electrically-induced phase difference Δθ_(e) (ω_(e)) between voltage (V) and current (I) within the energy device from applying the electrical excitations; measuring a mechanically-induced phase difference Δθ_(e) (ω_(m)) between voltage (V) and current (I) within the energy device from applying the mechanical excitations; and deducing the state-of-health of the energy device by comparing the electrically-induced phase difference Δθ_(e) (ω_(e)) with the mechanically-induced phase difference Δθ_(e) (ω_(m)).

This same method can be used to provide accurate run-time and electrical performance information and life time predictions for the energy device. In one embodiment of the invention, one applies this testing method and determines a baseline phase difference between current and voltage. Then, at a later time, when it is desired to test the health of the energy device, one again applies this testing method and compares the later phase difference with the baseline phase difference. The later phase difference diminishing from the baseline phase difference indicates energy device aging, and a decline in the energy device's electrochemical health, roughly linearly with the diminution. Using various prediction algorithms' and prediction techniques this method and device can be used to accurately predict the energy devices residual energy capacity available for use as well as the devices remaining life time.

BRIEF DESCRIPTION OF THE DRAWINGS

The features of the invention believed to be novel are set forth in the appended claims. The invention, however, together with further objects and advantages thereof, may best be understood by reference to the following description taken in conjunction with the accompanying drawing(s) summarized below.

FIG. 1 is a graph illustrating the phase relationship between voltage (V) and current (I) in an ideal capacitor.

FIG. 2 illustrates a model of a flooded type lead-acid battery.

FIG. 3 is a graph illustrating the phase relationship between voltage (V) and current (I) using a single-frequency EIS technique.

FIG. 4 is a graph illustrating the phase relationship between voltage (V) and current (I) using a dual-frequency EIS technique.

FIG. 5 is a graph illustrating the phase relationship between voltage (V) and current (I) using an ultrasonic technique.

FIG. 6 is a graph illustrating in general terms, true state of health measurement in accordance with various embodiments of the invention based on voltage (V) versus current (I) phase differences using electrical and mechanical excitations.

DETAILED DESCRIPTION

This disclosure will elucidate the novelty, inventive step and industrial viability of this method and device by illustrating the drawbacks of current SoH measurement techniques on lead acid battery chemistry, however, it should be understood that this technique can be applied to any electrochemical energy storage and/or conversion device; including, certain types of fuel cells, solar cells and capacitors (i.e. hybrid and super capacitor technologies).

In general, many of the inherent limitations of typical electrochemical energy storage and/or conversion devices are the result of kinetic deficiencies, i.e. low diffusion and migration rates of the chemical reaction product ions (i.e., ionic species). In fact, these same kinetic deficiencies are the main differentiator between the performance of an ideal capacitor and an electrochemical energy storage device (e.g., battery). Certain types of fuel cells, such as: Polymer Electrolyte Membrane (PEM) also have ionic kinetic deficiencies for which this method and device can be exploited to gain additional electrochemical health information which will help provide accurate SoH measurement and prediction methods. In an ideal capacitor there is no chemical reaction taking place and therefore the ionic species kinetic transport limitations are non-existent and do not impact this storage device's performance, i.e. almost instantaneous charge time and almost infinite cycle life. However, ideal capacitors do not have the energy storage capability that would be useful for large critical energy storage back-up applications. In fact, their storage capacities are orders of magnitudes less than that of even the most limited electrochemical storage device (e.g., battery). There are other types of capacitors on the market and in development in which there are chemical reactions taking place and which therefore suffer from ionic kinetic limitations, however; in this disclosure we will utilize the behavior of this ideal capacitor as an example, without implied limitation, to elucidate the effectiveness of the disclosed novel and innovative EIS SoH monitoring method and device.

The Phase Relationship Between Voltage (V) and Current (I) Across an Ideal Capacitor

FIG. 1 is a phase diagram that will be used throughout the remainder of this disclosure to provide insight into the novelty and inventiveness of this method and device. In this disclosure the use of the term ‘ideal’ in reference to the ideal capacitor is used to elucidate the lack of ionic kinetic deficiencies that plague many other types of electrochemical energy storage and/or energy conversion devices. In this diagram the phase relationship between Voltage (V) and Current (I) when applied to an ideal capacitor is depicted. A very small amplitude sinusoidal voltage is applied across the terminals of the capacitor and the resultant phase difference (Δθ) of the current response is measured. In the case of an ideal capacitor this Δθ difference is exactly 90°. One may ask that if we applied V first then why does the current (I) show up first in the phase diagram. This is not intuitive, so it will be discussed briefly. When the electric potential is felt across the capacitor electrodes, electrons start to accumulate on the negative electrode and then they essentially push electrons off the positive electrode, thus leaving that electrode positively charged, and this is when the V across the capacitor electrodes begins to grow. So, when all the electron flow stops (point A in FIG. 1), the maximum V is measured across the capacitor terminals.

The phase relationship depicted in FIG. 1 is that of an ideal capacitor. The only charge carriers involved are electrons and there are no chemical reactions happening within that capacitor and therefore there are no ionic kinetic limitations or deficiencies. Clearly, if this same V signal was applied across the terminals of a battery, an electrochemical device which does undergo chemical reaction steps, the phase difference Δθ between V and I will be very different. This new phase difference will be a result of many factors, but, the most significant difference will be due to the introduction of a new type of charge carrier, the chemical product ion charge carriers. These ion charge carriers are much slower than electrons and their rates are dependent on many different factors which will be discussed below.

The Single Electrical Excitation Frequency (ω_(e)) Technique

Early attempts at determining the health of a battery without invoking disruptive and costly test discharges fell short. With this technique, a small sinusoidal electrical excitation voltage is applied across the battery's terminal and the resultant phase difference Δθ of the response current is measured. The clear advantage of this method is that it can be used without taking the critical battery back-up off-line and it is a very fast test. If this method of SoH determination was accurate then it would obviate the need to take the battery pack out of service, and it would eliminate the resultant operational restrictions discussed above. The phase difference Δθ is correlated to the battery's internal impedance (Z_(INT)), thus it is an indirect measure of the battery's opposition to the flow of energy and this can be related to the health of the battery. In order to examine this technique closer a typical and accepted model of a flooded type lead-acid battery which will be used and it is depicted in FIG. 2, where:

-   -   R_(m)=Metallic Resistance (post, alloy, grid, and gridpaste)     -   R_(el)=Electrochemical Resistance (paste, seperator, and         electrolyte)

In FIG. 2, the top depiction models the physical components that offer an opposition to the flow of energy through the battery. The equivalent circuit (the bottom depiction) has grouped these components in order to better describe the battery's behavior as it relates to different types of charge carriers (i.e. ions and electrons). The capacitance symbol (C_(b)) represents the storage capacity of the battery. The metallic resistance (R_(m)) represents the opposition to the flow to electron charge carriers, and the electrochemical resistance (R_(el)) represents the opposition to the flow of ion charge carriers. This breakout of charge carrier oppositions will aid in the following description:

-   -   (1) Complex Impedance of Battery (Z_(INT))

$Z_{INT} = {\left( {{R_{m}\left( \omega_{e} \right)} + \frac{\frac{1}{R_{el}\left( \omega_{e} \right)}}{\frac{1}{R_{el}^{2}\left( \omega_{e} \right)} + \left( {\omega_{e}{C_{b}\left( \omega_{e} \right)}} \right)^{2}}} \right) - {j\left( \frac{\omega_{e}{C_{b}\left( \omega_{e} \right)}}{\frac{1}{R_{el}^{2}} + \left( {\omega_{e}{C_{b}\left( \omega_{e} \right)}} \right)^{2}} \right)}}$

where;

-   -   C_(b)=Battery Capacitance     -   ω_(e)=electrical excitation Frequency

Equation 1 above was derived by performing a simple AC analysis of the equivalent circuit in FIG. 2. Z_(int) is an expression for the total complex impedance of the battery, which represents the opposition to the flow of energy to an AC electrical excitation signal with a frequency (ω_(e)). As a note: If a DC signal is applied, then ω_(e)=0, and the Z_(int) expression simplifies to Z_(int)=R_(m)+R_(el).

-   -   (2) Voltage (V) and Current (I) phasor relationship

$\frac{V\left. \sqrt{}\theta_{e} \right.}{I} = Z_{INT}$

Where;

-   -   θ_(e)=Phase Angle due to electrical frequency excitation

Equation 2 above is an expression of the phase relationship between V, I and Z_(INT). This is merely Ohm's Law rearranged and in its complex form. The resultant phase angle (θ_(e)) is presented with sub-script (e), and this is associated with the electrical excitation frequency (ω_(e)) presented in Equation (1). These subscripts are not typically shown, but this granularity will be useful in demonstrating the novelty of the method being disclosed.

As shown in Equations (2), the phase angle (θ_(e)) can be correlated to Z_(INT). It should be recognized that if Z_(INT) changes its response to an AC excitation signal, for any reason, then its resultant complex impedance response, at a certain frequency (ω_(e)), will change as well and will result in a different phase angle (θ_(e)). If it is assumed that the V, I and ω_(e) between various measurements are the same, then the resultant phase angle (θ_(e)), the measured parameter, will be dependent on Z_(INT) only. Thus, the phase angle response between V and I at a specific frequency (ω_(e)) can be correlated to the complex impedance Z_(INT) of a battery at the time the measurement is taken and in theory can be used to indirectly determine the health of a battery. However, in practice, this method is far from accurate as will be discussed later in this disclosure.

The Battery's Complex Impedance Z_(INT)

As the battery ages the available active material (energy capacity—Ahrs) decreases, lead sulfate and other tertiary compounds crystallize and form unwanted bonds, causing non-rechargeable masses to form on, around and within the electrode paste, the charge acceptance declines, the cyclic efficiency falls off, the power density declines, and eventually the battery cannot serve its intended function. Each of the components in FIG. 2 are effected by battery aging phenomenon in a different way and each of their contributions will impact the battery's complex impedance (Z_(INT)) in a different way. This will show up in a difference in the measured phase angle (θ_(e)) between the V and I at a specific frequency (ω_(e)). In Equation 1 above, it is clear that R_(m), R_(el), C_(b) and ω_(e) all affect the complex impedance response to an electrical excitation signal.

In addition to aging phenomenon, there other battery conditions that also affect the complex impedance response of the battery, such as; the State-of-Charge (SoC), the battery's internal Temperature (T), the battery design, and the battery manufacturing tolerances and process. These conditions affect the battery components, and their contribution to Z_(INT), and this also results in a difference in the measured phase angle (θ_(e)). A closer look into how each of these components are affected by these factors is necessary in order to better understand the problems with the accuracy of the single frequency excitation measuring technique. The key to understanding the limitations on the accuracy which is possible with the single frequency excitation method is to understand the variability that is introduced into this measurement technique by the different components in the battery model depicted in FIG. 2. Each of these components are impacted differently by the battery's age, SoC, temperature and its design (material and construction), manufacturing and process tolerances. In order for this technique to be an accurate and trusted SoH monitor it must reveal information about the true electrochemical health of the battery. Its measurement result cannot be altered by battery conditions such as SoC, temperature and its design or it will not supply the proper information to the battery user in order for them to make critical and costly decisions about their battery pack.

How does Battery Aging Affect Z_(INT)?

In general, R_(post), R_(alloy) and R_(grid), which are all components of R_(m), are not significantly impacted by battery aging. However, the following components are affected:

-   -   a) R_(Gridpaste), a component of R_(m). During the first several         cycles of a lead acid battery the bonding between the grid and         paste tends to get stronger, after that, the bonds tend to break         down. Thus, the contribution of this component to the resultant         phase angle (θ_(e)) will cause it to go one way at the beginning         of life and then go the other way as the battery ages, further         adding to the variability of this type of measurement technique.     -   b) R_(paste), a component of R_(el), increases as the battery         ages and this is mainly a result of unwanted tertiary compounds,         i.e. lead oxides and basic lead sulfates (tribasic and         tetrabasic), precipitating out of solution in unbalanced         alkalinity conditions during the discharge. These tertiary         compounds interact and establish bonding relationships with the         rechargeable PbSO₄ volume and the electrode, on the surface of         the electrode and within its pores. These interactions (i.e.         bonds) are mainly of the covalent or even metallic type, but         regardless of the type, they are undesired because they         negatively impact the storage device's capacity. This         effectively limits the flow and subsequent ionic transfer of the         electrolyte's chemical ion products and presents itself as an         increase to the opposition of the flow of energy, that is, an         increase in Z_(INT). This is essentially a decrease in the         active material utilization (measure of the battery's         capacity—Ahrs). One can see how this negatively impacts both the         discharge and charge of the battery. During discharge there is         less active material available to undergo energy conversion, and         during recharge any additional compounds and/or bonding that is         created during the discharge reaction must be undone during the         recharge reaction (i.e., decrease in charge acceptance). The         resultant contribution of this component tends to increase the         phase angle (θ_(e)) as the battery ages thus further adding to         the variability of this measurement technique.     -   c) R_(sep), a component of R_(el), tends to break down as the         battery ages and this also affects the resultant phase angle         (θ_(e)), thus further adding to the variability of this         measurement technique.     -   d) R_(electrolyte), a component of R_(el), is rather complicated         since there are so many different types of electrolytes         currently used in the industry: Liquid, Gel, and Absorbent Glass         Matting (AGM), etc. The list is large and the aging mechanism of         each is different, some worst then others; however, in general,         the contribution of this component to the resultant phase angle         (θ_(e)) will be affected by battery age thus further adding to         the variability of this type of measurement technique.     -   e) C_(b), a component of the battery model in FIG. 2 is also         affected by aging. Since C_(b) is a measure of the battery's         internal capacitance the same aging effects described for         R_(Paste) and R_(Electrolyte) also affect the battery's capacity         to supply energy. Generally the variability of the measurement         technique due to this component increases as the battery ages.         How does Battery State-of-Charge (SoC) Affect Z_(INT)?

In general, R_(post), R_(alloy), R_(grid), R_(gridpaste), which are all components of R_(m), and R_(sep), components of R_(el), are not significantly impacted by the battery's SoC. The following components are affected:

-   -   a) R_(electrolyte), a component of R_(el), is highly dependent         on the battery's SoC. The opposition to the flow of energy, i.e.         ions, increases as the concentration of the SO₄ ⁻⁻ in the         electrolyte solution goes down. Effectively, there are less SO₄         ⁻⁻ to react with the active material ions (Pb⁺⁺) and this         presents itself as an increase to the opposition of flow.         Therefore, the contribution of this battery condition will         change the resultant phase angle (θ_(e)), thus, further adding         to the variability of this measurement technique. This         variability also increases as the battery discharge continues.     -   b) R_(paste), a component of R_(el), is dependent on the         battery's SoC. If the battery cell is not fully charged than         there is lead sulfate layer on the surface and within the pores         of the plate. Thus, this layer increases the opposition to the         flow of ions. Therefore, the contribution of this battery         condition will change the resultant phase angle (θ_(e)), thus         further adding to the variability of this measurement technique.         This component of R_(el) is also the component that causes the         battery's terminal voltage to go down when the battery reaches         very low SoC. During the first major portion of the discharge of         a lead acid battery the output voltage is relatively constant,         however, near the end of discharge the voltage rapidly         decreases. As Z_(INT) rapidly increases this marks the end of         discharge. Thus, the contribution of this battery condition will         also change the resultant phase angle (θ_(e)) and the         variability of this measurement technique will increase as the         discharge continues     -   c) C_(b), a component of the battery model in FIG. 2 is also         affected by SoC. If the battery cell is not fully charged than         there is lead sulfate layer on the surface and within the pores         of the plate. Thus, this layer alters the permittivity of the         cell thus affecting the representative capacitance of the cell.         Therefore, the contribution of this battery condition will         change the resultant phase angle (θ_(e)), thus further adding to         the variability of this measurement technique.         How does Battery Temperature Affect Z_(INT)?

In general, R_(Sep), a component of R_(el), could be temperature sensitive, but, in general it is not over the temperature ranges of interest, therefore in most practical situations this component will have minimal impact on the resultant phase angle (θ_(e)). The effects of the contribution of temperature on C_(b) is minimal compared to the effects noted below. The following components are affected by temperature:

-   -   a) R_(Grid), R_(Alloy), R_(post), and R_(Gridpaste), components         of R_(m), are temperature sensitive. Therefore, a majority of         the metallic opposition (R_(m)) to the flow of energy, i.e.         electron flow, depends of the battery's temperature when the         measurement is taken. Thus, this battery condition affects the         contribution that these components add to the resultant phase         angle (θ_(e)), further adding to the variability of this         measurement technique.     -   b) R_(electrolyte), and R_(Paste), components of R_(el), are         highly sensitive to temperature. Anytime one heats up a chemical         reaction the reaction rate will increase. Thus, this battery         condition affects the contribution that these components add to         the resultant phase angle (θ_(e)), further adding to the         variability of this measurement technique.         How do Battery Design, and/or Manufacturing Tolerances and         Process Affect Z_(INT)?

In general, all the components of the battery are affected by its design/manufacturing tolerances and process. However, some are not as significant as others.

-   -   a) The components of R_(m) are subject to variations based on         different battery designs and the materials used. More         importantly, these components are highly susceptible to         manufacturing tolerances and processes, i.e. various         welding/bonding techniques. Each manufacturer makes the battery         slightly different, therefore the complex impedance response for         the same battery chemistry will be different. Just because all         the batteries are of the same lot does not mean their electrical         response characteristics are the same. Manufacturing batteries         is a messy and highly variable process. Particularly, every         battery will have its own impedance signature, thus affecting         the resultant phase angle (θ_(e)) and further adding to the         variability of this measurement technique.     -   b) The paste, separator and electrolyte are tightly controlled         via industry regulations. R_(paste), R_(sep) and         R_(electrolyte), components of R_(el), and are less susceptible         to design, manufacturing tolerances and/or process variations.         Nonetheless, the complex impedance response and thus the         resultant phase angle (θ_(e)) are affected by this factor, thus,         adding to the variability of this measurement technique.     -   c) Clearly C_(b) is very dependent on the design and         manufacturing process and therefore this battery condition         affects the contribution that this component adds to the         resultant phase angle (θ_(e)), further adding to the variability         of this measurement technique.

It is worth noting that since each battery is a bit different, i.e. has a different electrical excitation frequency response, the actual EIS equipment that different manufacturers design also add to the variability of this measurement technique. So even a battery from the same lot and under tight industry regulation will respond differently and different conclusions will be drawn on the health of the battery, based on the resultant phase angle (θ_(e)), because the EIS manufacturer uses a different electrical excitation frequency. Manufacturers have not settled on the same (the right) electrical excitation frequency and this alone adds variability to this measurement technique.

In summary, the single frequency measurement technique comes with significant measurement variability and this takes away from its accuracy. Therefore, any results from this technique will remain in question and its ability to eliminate the costly, restrictive and time consuming test discharge has fallen short. At best, this technique may be useful in identifying imminent failure, i.e. a rogue cell, but this is too little, too late. Consequently, based on the inherent variability in this measurement technique, if it were the only indicator the utility used in determining the health of a battery pack then costly ‘premature’ battery pack replacement would be likely based on uncertainty alone.

FIG. 3 shows the resultant phase angle (θ_(e)) using the single excitation frequency EIS technique. The green shaded area is representative of the phase difference from that an ideal capacitor would present at the same frequency (recall from FIG. 1 that for an ideal capacitor V reaches a maximum or minimum just as I is zero, and vice versa). The difference is based on the battery's inherent kinetic transport limitations, i.e. slower ion flow, but, it is also comprised of the variability and uncertainly that was a result of the change in the battery's complex impedance response from the effects of aging, SoC, temperature, and battery design, manufacturing tolerances and/or processes. Therefore, it should be clear that only part of the green shaded area truly represents the battery's SoH, it's age, and the rest represents inherent features of the battery even for a brand new battery with excellent SoH.

One could surmise that each of the EIS equipment suppliers could determine the complex impedance response of each and every battery that is available and use this beginning of life measurement to yield better results. This equipment manufacturer could keep an extensive database so when their equipment was used it would give a better indication of that battery's health. However, this is a daunting task to say the least; first battery manufacturers are not that forthcoming with battery data and they would most likely not provide this initial measurement, secondly there are way too many types of batteries available to keep track of, and finally this would only remove some of the variability's discussed above. Therefore, this attempt would significantly complicate and add cost to this measurement technique and the benefits for undertaking this approach would be obviated by the remaining variability (inaccuracy) that was still in the measurement.

True State-of Health (SoH) Meaning

SoH has many different meanings and is often misinterpreted. True SoH measurement is supposed to give the user a better understanding of the chemical potential energy (Joules) that is available for conversion to usable electrical energy (Ahrs) when the battery is called upon. If the measurement technique used is influenced by factors such as SoC, temperature, and battery design, manufacturing tolerances and/or process, as is the single frequency excitation technique, then the user will not be presented with the right information to make a decision whether or not to replace the battery pack. This uncertainty is costly to say the least and thus this is why utility companies still to this day rely on costly, restrictive, and time consuming periodic test discharges.

What we really want to know from a SoH measuring technique is the battery pack's true electrochemical health. The SoC, temperature and battery design, manufacturing tolerances and/or process, and their contributions to the battery's complex impedance response, and resultant phase angle response, does not tell us whether or not our battery pack is healthy and we need a way to factor these out. What we really need to know in order to make an informed decision about our battery's health is how old is it, the battery's age in electrical terms! Clearly, we want to know how our battery pack has aged based on the regimen we have put it through. So, in order for a measuring technique to be informative we need to eliminate the variability that is introduced by its components based on SoC, temperature, and battery design, manufacturing tolerances and processes and we also need to eliminate the variability introduce by different EIS measuring equipment. As discussed previously, each EIS device manufacturer uses different electrical excitation frequencies and this alone introduces measurement variability. If one could get all the manufactures to use the same electrical excitation frequencies then this variability could be substantially reduced, but this is highly unlikely. The resultant measurement needs to tell us about how the age of the battery is going to impact its ability to serve its intended function when called upon, Joules to AHrs conversion. Thus, the measuring technique needs to somehow eliminate the variability's discussed above.

From FIG. 3, the ultimate goal is to reduce the green shaded area, which would mean that some of the variability introduced by factors other than battery aging would have to be eliminated and the resultant measurement would be more indicative of the battery's true SoH, i.e. its electrochemical health.

Dual Electrical Excitation Frequency (ω_(e) ¹ and ω_(e) ²) Attempt

As a natural response to the shortcomings of the single frequency measurement method, an attempt was made to use two electrical excitation frequencies, ω_(e) ¹ and ω_(e) ² respectively. The goal here was to extract additional electrochemical health information, by eliminating some of the variations introduced by non-relevant factors, by exciting the battery with two different electrical excitation frequencies, close in time. However, this method also fell short. The following relationships illustrate this concept:

-   -   (3) Complex Impedance of Battery (Z_(INT)), with ω_(e) ¹ and         ω_(e) ²

$Z_{INT}^{1} = {\left( {{R_{m}\left( \omega_{e}^{1} \right)} + \frac{\frac{1}{R_{el}\left( \omega_{e}^{1} \right)}}{\frac{1}{R_{el}^{2}\left( \omega_{e}^{1} \right)} + \left( {\omega_{e}^{1}{C_{b}\left( \omega_{e}^{1} \right)}} \right)^{2}}} \right) - {j\left( \frac{\omega_{e}^{1}{C_{b}\left( \omega_{e}^{1} \right)}}{\frac{1}{R_{el}^{2}\left( \omega_{e}^{1} \right)} + \left( {\omega_{e}^{1}{C_{b}\left( \omega_{e}^{1} \right)}} \right)^{2}} \right)}}$ $Z_{INT}^{2} = {\left( {{R_{m}\left( \omega_{e}^{2} \right)} + \frac{\frac{1}{R_{el}\left( \omega_{e}^{2} \right)}}{\frac{1}{R_{el}^{2}\left( \omega_{e}^{2} \right)} + \left( {\omega_{e}^{2}{C_{b}\left( \omega_{e}^{2} \right)}} \right)^{2}}} \right) - {j\left( \frac{\omega_{e}^{2}{C_{b}\left( \omega_{e}^{2} \right)}}{\frac{1}{R_{el}^{2}} + \left( {\omega_{e}^{2}{C_{b}\left( \omega_{e}^{2} \right)}} \right)^{2}} \right)}}$

Where;

-   -   ω_(e) ¹=1^(st) electrical excitation Frequency     -   ω_(e) ²=2^(nd) electrical excitation Frequency     -   (4) Voltage (V) and Current (I) phasor relationship, with ω_(e)         ¹ and ω_(e) ²

$\frac{V\left. \sqrt{}\theta_{e}^{1} \right.}{I} = Z_{INT}^{1}$ $\frac{V\left. \sqrt{}\theta_{e}^{2} \right.}{I} = Z_{INT}^{2}$

Where;

-   -   θ_(e) ¹=Phase Angle difference between V and I due to 1^(st)         electrical frequency excitation     -   θ_(e) ²=Phase Angle difference between V and I due to 2^(nd)         electrical frequency excitation

In Equation 4, there are two different complex internal impedances responses, Z_(INT) ¹ and Z_(INT) ². The battery responds differently because it is being excited by a different electrical excitation frequency and the resultant complex internal impedance response to each is different. Some of the components that make up the battery's complex internal impedance (Z_(INT)) present a different opposition to the various frequencies with which they are excited. The components that make up R_(m) (the metallic opposition, electrons) are not very frequency sensitive at the frequencies of interest. However, the components that make up R_(el) (the electrochemical opposition, ions) are sensitive at all frequencies. Additionally, the capacitive component, in FIG. 1, and its contribution to the complex impedance response of the battery is also frequency sensitive. Thus, it should be clear, that the two expressions for Z_(INT) in Equation 4, Z¹ _(INT) and Z² _(INT), will have a different response to different excitation frequencies, ω₁ and ω₂ respectively. But, recall that the main problem with the single frequency method was that the complex impedance response, and it's resultant phase angle, was highly dependent on the age, SoC, temperature and various design, manufacturing tolerances and/or process. So, in order for this approach to be effective, and to instill confidence in this electronic battery SoH measurement technique, it would have to eliminate all of these dependencies so that the resultant measurement is truly indicative of the cell's age.

This measurement technique takes the difference of the resultant phase angle response, θ_(e) ¹ and θ_(e) ², measured at the terminals of the battery, and this provides additional insight in to the battery's internals and this aids in determining the true SoH of the battery. One would think that this two frequency approach is far superior to the single frequency approach; however, it is the intent of this disclosure to show otherwise. Just using two electrical excitation frequencies taken close in time is not enough to make this method of measurement acceptable in the industry. There are other variations of this concept that have been explored, e.g. scanning through a wide range of frequencies and measuring several resultant responses, all of the electrical excitation type. These attempts have merit but they are also limited and they too have several variations in the response which make these measurements unreliable as well. The added scanning significantly complicates the required electronics and greatly raises the cost. One must keep in mind that the utility or end user is not going to pay more for a battery monitoring system then the battery itself. This would not make much business sense; why not just replace the battery when it is suspect. So, the goal is to eliminate the variation without complicating the required electronics and increasing the cost.

-   -   (5) The battery internal impedance difference ΔZ_(INT) based         on (3) is:

${\Delta \; Z_{INT}^{e}} = {{Z_{INT}^{e\; 1} - Z_{INT}^{e\; 2}} = {\left( {\frac{\frac{1}{R_{el}\left( \omega_{e}^{1} \right)}}{\frac{1}{R_{el}^{2}\left( \omega_{e}^{1} \right)} + \left( {\omega_{e}^{1}{C_{b}\left( \omega_{e}^{1} \right)}} \right)^{2}} - \frac{\frac{1}{R_{el}\left( \omega_{e}^{2} \right)}}{\frac{1}{R_{el}^{2}\left( \omega_{e}^{2} \right)} + \left( {\omega_{e}^{2}{C_{b}\left( \omega_{e}^{2} \right)}} \right)^{2}}} \right) - {{j\left( {\frac{\omega_{e}^{1}{C_{b}\left( \omega_{e}^{1} \right)}}{\frac{1}{R_{el}^{2}\left( \omega_{e}^{1} \right)} + \left( {\omega_{e}^{1}{C_{b}\left( \omega_{e}^{1} \right)}} \right)^{2}} - \frac{\omega_{e}^{2}{C_{b}\left( \omega_{e}^{2} \right)}}{\frac{1}{R_{el}^{2}\left( \omega_{e}^{2} \right)} + \left( {\omega_{e}^{2}{C_{b}\left( \omega_{e}^{2} \right)}} \right)^{2}}} \right)}.}}}$

The impedance difference in Equation 5, although tedious, represents the advantages of the dual electrical frequency approach. Most importantly, this complex impedance response does not contain R_(m) which drops out when taking this difference, and as will be recalled from the single frequency method description this term was present in the complex impedance (Z_(int)) and added variability to the measured phase angle (θ_(e)). Only one component of R_(m), R_(gridpaste), was affected by age, but all the components of R_(m); R_(post), R_(alloy), R_(grid) and R_(gridpaste) were all affected by temperature and the battery design, manufacturing tolerance and/or process. Thus, eliminating this term from the complex impedance expression was a step in the right direction. That is, this dual frequency excitation difference method eliminated the effects that two of the battery's conditions, i.e. temperature and the battery, its design, manufacturing tolerance and/or processes, had on R_(m), and thus reduced this contribution to the variability of the measurement.

By exciting the battery close in time and taking the difference of the phase response of two distinct electrical excitation frequencies two of the large contributors to the variation of the measurement have been eliminated. Hence, the variability in this measurement technique has been reduced, and closer approaches a measurement that is more representative of the true SoH of the battery, its age. Essentially, since the temperature did not change and the battery design is the same between the two measurements their contribution to the variability in the measurement is eliminated. FIG. 4 shows that the green shaded area in FIG. 3 above has been reduced by this dual frequency approach. The red shaded area represents the variability that was removed by eliminating the effects of temperature and battery design, manufacturing tolerances and/or process on R_(m).

However, from Equation 5 it is shown that Z_(INT) is still dependent on R_(el), C_(b) and ω_(e). As described earlier, R_(el) is affected by the battery's age, SoC, temperature and battery design, and manufacturing tolerances and/or processes. So even with this dual frequency excitation difference approach there will still be significant variability in the measurement and it is not all due to the battery's age. The main goal of the dual frequency approach was to eliminate the variation in measurement so that it better represented the true health of the battery, however, in doing so, this measurement technique actually added variability into the measurement results. Since, the components that make up R_(el) will behave differently to different electrical excitation frequencies this behavior manifests itself the difference measurement. Thus, new unknowns have been added to the measurement technique that were non-existent in the single electrical excitation technique. Therefore, the selection of the second excitation frequency alone will yield different results in the measured phase angle difference, Δθ_(e).

It is clear that even though the dual frequency approach helped, and was superior to the single frequency approach, it still has limitations in its ability to accurately measure the SoH of a battery. The variations induced by different battery types and material, i.e. paste, separator and electrolyte, aging, SoC, temperature, are still present within the ΔZ_(INT) expression in Equation 5 and they are a direct result of the variations in R_(el) and its frequency response to two distinct electrical excitation frequencies.

It is worthy to note that much study and exploration has been expended on the electrical frequency selection for this type of measurement method. From Equation 5 it is evident that the value of both ω_(e) ¹ and ω_(e) ² will affect the phase difference, Δθ_(e). Meaning, if two identical batteries are tested and all other parameters remained equal one would get two different Δθ_(e), if one used two differently supplied EIS equipment. This is not visible to the end-user because the equipment manufactures draw their own conclusions on what the actual SoH is based on the phase difference response seen by their system. The real goal was to reduce variability in the measurement of SoH, in order to get the true electrochemical health of the battery, now we have even furthered the level of variation: the EIS equipment supplier's selection of two electrical excitation frequencies. So, we eliminated some of the variability of the single excitation technique but added some variability based on our enhanced measurement technique. Thus, this tends to keep the green shaded area in FIG. 4 larger than it should be. Consequently, it is not surprising that one EIS vendor may say that their system is superior based on frequency selection, but, as we have disclosed, all EIS equipment using this electrical measurement method is subject to the variations brought on by R_(el). Albeit, the variation may be less from one frequency selection to the next, however, this difference is marginal at best. Certainly not a ‘real’ performance differentiator in the market place!

The Mechanical Excitation Frequency (ω_(m)) and Single Electrical Excitation Frequency (ω_(e)) Method

As shown the single and dual electrical excitation frequency approaches both have their advantages and disadvantages. The evolution of EIS techniques/technologies has grown in complexity and cost but the performance, i.e. their ability to accurately measure SoH, has not grown in proportion. It is disclosed here that there is a much more accurate method for determining the actual SoH of a battery. The principle behind this method is to introduce a mechanical excitation frequency ω_(m) as is disclosed in U.S. Pat. No. 7,592,094. This prior art demonstrates that the introduction of a mechanical excitation frequency ω_(m) to an energy storage device makes the device operate much more efficiently and will substantially minimize the storage device's degradation. This disclosure is focused on the advantages that mechanical excitation can additionally have on one's ability to accurately determine the battery's true SoH, its age. This adds value to the use of mechanical excitation in the energy storage market place, meaning, if the utility is taking advantage of the performance and cycle life improvements that mechanical excitation offers then it makes commercial business sense to take advantage of this novel and industrious method of SoH determination. The combination of both mechanical excitation for improved performance and life and a better and more accurate SoH measurement technique will change the face of energy storage as we currently know it.

It will be demonstrated that the variations brought on by R_(el), in the single and dual electrical excitation frequency can be all but eliminated by the proper application of a mechanical excitation frequency (ω_(m)). As discussed above, the expression for ΔZ_(INT) in Equation 5 was effected by the battery's SoC, temperature and various design, manufacturing tolerances and/or process as well as the device's age, i.e. SoH. When testing a string of batteries with the same electrical excitation measurement device, the user will get a different Δθ_(e), i.e. SoH, reading for each battery in the string. The goal of any measurement technique is to ensure that the measured Δθ_(e) is indicative of the battery's true SoH, and of that alone.

To illustrate, take a string of batteries and envision that we are going to test each one individually. Each battery in the string will be at a different age, even if the string was comprised of batteries from the same lot and placed into service at the same time. The problem is that each battery cell charges and discharges at different rates throughout its life and this is mainly due to differences in its internal impedance (Z_(INT)) and imbalances between cells in the string. These are the same differences that make the current electronic excitation means of determining the batteries true SoH very inaccurate and unreliable. Therefore, each cell will be at a different cycle life, a.k.a. age. Often, one confuses cycle life with calendar life but the differentiation is critical.

Since it was already determined that R_(el)'s magnitude is sensitive to aging mechanisms, part of each cell's complex impedance response to an electrical excitation frequency, and the corresponding Δθ_(e), is due to age. The question is how much of the measured phase difference can be attributed to age? After all, this is what everyone is looking for since the battery's true SoH is, or more correctly, should be, a measure of its electrochemical health, e.g. its ability to convert chemical potential energy to electrical energy then back again. So, what we really want to know is how the paste, separator and electrolyte will interact, i.e. the energy conversion process, when called upon for service. The effects from the battery's SoC, temperature, and design, manufacturing tolerance and/or process should not be present in the measurement. Even though these parameters play a part in what the battery will supply, i.e. capacity—Ahrs, for its given state, one should not mislead one into prematurely replacing a battery pack based on the variability they add to the chosen measurement technique. All we really want in the measurement result are the effects that are due to the battery's age! R_(el)'s magnitude in response to an electrical excitation frequency, or the difference in the response from two electrical excitation frequencies, alone is descriptive, but it is also very deceiving and this inherent uncertainty directly effects a utility's bottom line, i.e. cost! Clearly, if a battery pack's SoH is in question based on the results of the applied measurement technique then this will drive the utility to performing a costly and operational restrictive test discharge which, as explained earlier, actually detracts from the battery pack's useful life. Therefore, the purpose of this novel, inventive and industrious mechanical excitation technique is to remove this inherent uncertainty in the electronic measurement technique, thereby, obviating the need to conduct frequent and costly test discharges.

R_(el) is comprised of three components; R_(Paste), R_(Sep) and R_(electrolyte), their individual magnitude response to an electrical excitation frequency will be analyzed below:

-   -   1. R_(Paste), R_(Sep), R_(electrolyte) and C_(b) are sensitive         to aging, thus adding variability to resultant Δθ_(e)     -   2. R_(paste) and R_(electrolyte) are sensitive to temperature,         thus adding variability to resultant Δθ_(e)         -   a. As discussed earlier R_(Sep) could be sensitive to             temperature, but in general, in most practical             circumstances, it is not in the temperature range of             interest.     -   3. R_(electrolyte), R_(Paste) and C_(b) are sensitive to SoC,         thus adding variability to resultant Δθ_(e)     -   4. R_(Paste), R_(Sep), R_(electrolyte) and C_(b) are susceptible         to battery material and manufacturing tolerances and/or process,         thus adding variability to resultant Δθ_(e)

Thus as can be seen each of these component's magnitude response to an electrical excitation frequency are impacted by the conditions noted above, thus adding to the variability of the resultant measured phase angle θ_(e). Therefore the true effects of aging on the electrochemical health of the battery, i.e. SoH, are masked by the effects of temperature, SoC and the battery design, manufacturing tolerances and/or processes. It is also subject to the EIS equipment design and tolerances as well as this equipment supplier's interpretation of the electrical excitation measured response.

By measuring the internal impedance with and without a mechanical excitation frequency (ω_(m)) and only using a single electrical excitation frequency (ω_(e)) the influences of the temperature, SoC and battery design, manufacturing tolerances and/or processes on the measured phase angle θ_(e) can be eliminated and the resultant measured θ_(e) will be indicative of the true electrochemical health, i.e. SoH (age), of the battery. This mechanical measurement technique also reduces the response's susceptibility to the EIS equipment design and tolerance and substantially reduces variance introduced by the designer's interpretation of the result measured response.

-   -   (6) Battery Internal impedance difference, with ω_(e) and ω_(m),         ΔZ_(INT) ^(ω) ^(e) (ω_(m)), compare to equation (5) in which         both applied frequencies are electrical.

${{\Delta \; {Z_{INT}^{\omega_{e}}\left( \omega_{m} \right)}} = {{Z_{INT}^{\omega_{e}} - {Z_{INT}^{\omega_{e}}\left( \omega_{m} \right)}} = {\left( {\frac{\frac{1}{R_{el}\left( \omega_{e} \right)}}{\frac{1}{R_{el}^{2}\left( \omega_{e} \right)} + \left( {\omega_{e}{C_{b}\left( \omega_{e} \right)}} \right)^{2}} - \frac{\frac{1}{R_{el}\left( \omega_{m} \right)}}{\frac{1}{R_{el}^{2}\left( \omega_{m} \right)} + \left( {\omega_{e}{C_{b}\left( \omega_{m} \right)}} \right)^{2}}} \right) - {j\left( {\frac{\omega_{e}{C_{b}\left( \omega_{e} \right)}}{\frac{1}{R_{el}^{2}\left( \omega_{e} \right)} + \left( {\omega_{e}{C_{b}\left( \omega_{e} \right)}} \right)^{2}} - \frac{\omega_{e}{C_{b}\left( \omega_{m} \right)}}{\frac{1}{R_{el}^{2}\left( \omega_{m} \right)} + \left( {\omega_{e}{C_{b}\left( \omega_{m} \right)}} \right)^{2}}} \right)}}}},$

Constraints on Electrical Excitation

-   -   (7)

$\left( {\omega_{e}C_{b}} \right)^{2}\operatorname{>>}\frac{1}{R_{el}^{2}}$ ${{and}\left( {\omega_{e}{C_{b}\left( \omega_{m} \right)}} \right)}^{2}\operatorname{>>}\frac{1}{R_{el}^{2}\left( \omega_{m} \right)}$

Equation 6 can be manipulated and transformed by using constraints similar to that shown in Equation 7. One should keep in mind that there are other constraints on electrical excitation frequency selection that have been learned and applied by EIS equipment suppliers over decades of use and that these should also be considered during the implementation of this disclosed measurement technique. The mechanical excitation frequency selection is dependent on electrochemical constraints, i.e. chemical reactions and battery chemistry, which will be discussed later in this disclosure. The Z_(INT) expression in Equation 6 represents the battery's complex internal impedance and reveals its dependency on both the electrical and mechanical excitation frequency. Equation 8 and 9 below represent the polar form of the battery's internal impedance (Z_(INT)) and are used for illustrative purposes, and show that the resultant magnitude of Z_(INT) and the resultant phase shift θ_(e) between V and I is different when mechanical excitation ω_(m) is applied while taking a single electrical frequency measurement. The phase shift θ_(e) (Equation 8) is a result of electrical excitation and it is dependent on the mechanical excitation ω_(m) (Equation 9). The differences in the phase shift θ_(e) with and without mechanical excitation are due to the contributions of R_(el) and C_(b) in Equation 6 above.

Z_(INT) ^(ω) ^(e) ∠θ^(e) without mechanical excitation   (8)

Z_(INT) ^(ω) ^(e) (ω_(m))∠θ^(e)(ω_(m)) with mechanical excitation   (9)

This difference expression for Z_(INT), with and without mechanical excitation, reveals very pertinent electrochemical health information about the battery. From FIG. 5 the difference from the single electrical excitation phase response, and the single electrical excitation phase response with mechanical frequency excitation (the red shaded area), can be used to better ascertain the true SoH of the battery (the green shaded area). Essentially, this phase difference is with and without mechanical excitation is measured, and the resultant difference (green shaded area) is indicative of the electrochemical health of the battery.

The green shaded area (above the horizontal axis) in FIG. 5 does not have the variability that plagued the single electrical excitation method with the exception of that induced by the age of the device. These measurements are taken very closely in time in an environment in which the temperature is held substantially constant, therefore, temperature does not affect the change in the magnitude response of R_(el) with and without mechanical excitation, SoC does not affect the change in the magnitude response of R_(el) with and without mechanical excitation, and, finally the battery design, manufacturing tolerances and/or processes do not affect the change in magnitude response of R_(el) with and without mechanical excitation. These variability's are represented by the red shaded area (below the horizontal axis) in FIG. 5. All that affects the magnitude of R_(el) in this measurement technique is the effect of the mechanical excitation signal (ω_(m)) on the opposition to the flow of energy. Therefore, the battery's internal complex impedance response is not only dependent on the electrical excitation frequency (ω_(e)) it is also a function of the mechanical excitation frequency (ω_(m)).

By applying mechanical energy at the proper frequency and amplitude the kinetic transport deficiencies of the reaction ions are mitigated. The result is a reduction in R_(el) and this in turn will impact the complex impedance response and corresponding phase angle. This reduction in R_(el) is at the electrode-electrolyte interface and it is mainly due to the R_(paste) and R_(electrolyte) contribution. This is accomplished by the interaction of the mechanical wave and the subsequent altering of limiting pore sizes that are formed within the insulator layer (PbSO₄) that is found all over the electrode reaction area during discharge. The mechanical waveform is at a frequency which ensures that the pore dimensions are kept large enough to allow/ensure passage of the larger and less energetic reaction product ions (HSO₄ ⁻ and SO₄ ⁻⁻) so that their concentration levels and corresponding ionic fluxes are generally higher and they are more able to rapidly change, at distances that are closer to the electrode surface, in response to changes in the Pb⁺⁺ ion flux (which are based on electrode potential changes). If the mechanical energy does not reduce R_(el) then the electrochemical health of the cell is in question. As the battery ages the effects of mechanical energy on the magnitude of R_(el) will slowly diminish overtime, and will result in less of an impact on the battery's complex impedance response and therefore less of a change in the resultant phase angle.

It should also be seen that the variance introduced by the EIS equipment itself has been eliminated from the resultant complex impedance response and thus does not contribute to the variability in this measurement technique, as long as the electrical excitation frequency selection is sufficient so that the constraints illustrated in Equation 7 hold true.

The only thing that will influence the change in the magnitude response of R_(el) is the predetermined mechanical excitation signal. In this measurement technique the only variability introduced will be due to the effects of aging on R_(el) and C_(b). This is the goal of any EIS measurement technique.

True State-of-Health (SoH) Measurement using Single Electrical (ω_(e)) and Mechanical (ω_(m)) Frequency Excitation Technique

Since the variability in the single electrical frequency technique was eliminated by applying a mechanical excitation frequency, a simple and straightforward measurement of the true SoH of the battery is now possible. However, it should be understood that many different algorithms' and prediction models can be employed based on this disclosure. The underlying premise that makes these other modeling algorithms' and prediction methods deployable is that the application of mechanical excitation energy has provided additional information on the true electrochemical health of the electrochemical energy storage and/or energy conversion device. This measure is indicative of the battery's true electrochemical health and this will give the battery user reliable and unquestionable information about the battery's health, its age. The only change in the phase angle Δθ_(e) returned while the battery is under the influence of the proper mechanical excitation energy will be due to the battery components that are susceptible to aging phenomenon. As the battery ages these battery components will be less and less susceptible to the effects of the applied mechanical energy and the electrical phase change based on this mechanical excitation will become less and less prominent, i.e., Δθ_(e)(ω_(e))−Δθ_(e)(ω_(m))→0. Essentially, the green shaded area (above the horizontal axis) in FIG. 5 will get larger and larger as the battery undergoes aging phenomenon, with Δθ_(e)(ω_(m))→Δθ_(e)(ω_(e)). The magnitude of this phase angle change Δθ_(e)(ω_(m)), when mechanical excitation (ω_(m)) is applied, is dependent only on the contributions of the battery components, R_(Paste), R_(Sep), R_(electrolyte), and C_(b). As these components age, the mechanical energy (ω_(m)) will not have such a pronounced effect on their contributions to this electrical phase change Δθ_(e)(ω_(m)), and this will be indication that the electrochemical health of the battery is declining. A measurement of the Δθ_(e)(ω_(m)), while mechanical excitation (ω_(m)) is applied, will be the basis that all future measurements will be compared to. If at some future time in the battery's life, the measured Δθ_(e)(ω_(m)), while mechanical excitation (ω_(m)) is applied, is less than that of the previous measurement, then this difference is due to battery aging only. Eventually, when the battery is heavily aged, or near end of life, the mechanical excitation energy will not change the electrical phase angle θ_(e), i.e., the difference between these two phase differences Δθ_(e)(ω_(e))−Δθ_(e)(ω_(m))→0, and the battery pack, and/or cell, should be replaced.

Since the aging effect of the battery components in question will gradually reduce their respective contributions to the battery's internal complex impedance response, and correspondingly change the measured phase angle θ_(e), one can take measurements at predetermined times, chart the progression in the phase angle change, and actually predict when the battery, and/or cell, will reach the end of its useful life. In fact, if one knows the actual SoH of the device under test one can apply a SoH factor to a known SoC and be able to determine the available capacity of that energy storage device and/or the conversion efficiency of an energy conversion device. It should also be understood that the variability introduced by the majority of the battery's components in the single electrical frequency technique can also still yield useful information on the overall health of the battery.

Several different and useful variations of this prediction model can be incorporated, and this will then aid the battery user in making decisions about their battery pack's health and when imminent replacement is necessary. Several phase angle θ_(e)(ω_(m)) measurements, while mechanical excitation (ω_(m)) is applied, are taken and a projection of remaining battery life is determined. In one non-limiting embodiment of this disclosure, the age related contributions are projected using a straight-line approximation method and a prediction on remaining battery life is determined. Using various predictions algorithms and prediction techniques this method and device can also be used to accurately predict the device's residual energy capacity available for use as well as the device's remaining life time. While only the SoH determination method is described in this disclosure, it should be understood that the underlying premise of the measurement and prediction technique disclosed here is based on the availability of additional electrochemical health information that is made possible with the use of mechanical excitation energy with the current, and future, electrical excitation only measurement techniques, whether it be single frequency, dual, scanning, etc.

As a non-limiting example, FIG. 6 is the straight-line approximation method for a typical battery. Let us say for illustrative, non-limiting numbers, that when a given battery is new, a battery is tested, and as a baseline, it is established that Δθ≡Δθ_(e)(ω_(e))−Δθ_(e)(ω_(m))=80°. Over time as the battery ages, this 80° will diminish, which is to say, Δθ_(e)−Δθ_(e)(ω_(m)) will trend downward toward zero. As Δθ≡Δθ_(e)(ω_(e))−Δθ_(e)(ω_(m))→0°, this means that the battery is approaching it end of life. While the actual state of health may not be perfectly linear as illustrated in FIG. 6, i.e., while a Δθ_(e)(ω_(e))−Δθ_(e)(ω_(m))=40° reading which is halfway between 0° end of life and the 80° baseline may not mean that the battery is at exactly 50% of its useful life (i.e., while the straight, negatively-sloped line in FIG. 6 may not be exactly a line), one may use the linear scale in FIG. 6 as a good first order approximation of remaining battery life. FIG. 6 is analogous, then, to the indicator in an automobile which shows roughly the proportion of gasoline which remains in the tank. However, Δθ=Δθ_(e)(ω_(e))−Δθ_(e)(ω_(m))→0° is in all cases the measurement that indicates actual end of battery life, so the magnitude of Δθ_(e)(ω_(e))−Δθ_(e)(ω_(m)) in relation to a baseline-time phase difference (such as the 80° example above) may be taken as a rough linear gauge of battery SoH.

The knowledge possessed by someone of ordinary skill in the art at the time of this disclosure is understood to be part and parcel of this disclosure and is implicitly incorporated by reference herein, even if in the interest of economy express statements about the specific knowledge understood to be possessed by someone of ordinary skill are omitted from this disclosure. While reference may be made in this disclosure to the invention comprising a combination of a plurality of elements, it is also understood that this invention is regarded to comprise combinations which omit or exclude one or more of such elements, even if this omission or exclusion of an element or elements is not expressly stated herein, unless it is expressly stated herein that an element is essential to applicant's combination and cannot be omitted. It is further understood that the related prior art may include elements from which this invention may be distinguished by negative claim limitations, even without any express statement of such negative limitations herein. It is to be understood, between the positive statements of applicant's invention expressly stated herein, and the prior art and knowledge of the prior art by those of ordinary skill which is incorporated herein even if not expressly reproduced here for reasons of economy, that any and all such negative claim limitations supported by the prior art are also considered to be within the scope of this disclosure and its associated claims, even absent any express statement herein about any particular negative claim limitations.

Finally, while only certain preferred features of the invention have been illustrated and described, many modifications, changes and substitutions will occur to those skilled in the art. It is, therefore, to be understood that the appended claims are intended to cover all such modifications and changes as fall within the true spirit of the invention. 

I claim:
 1. A method for testing the state-of-health of an electrochemical energy device, comprising: applying electrical excitations to the energy device at a predetermined electrical excitation frequency ω_(e); applying mechanical excitations to the energy device at a predetermined mechanical excitation frequency ω_(m); measuring an electrically-induced phase difference Δθ_(e)(ω_(e)) between voltage (V) and current (I) within the energy device from applying said electrical excitations; measuring a mechanically-induced phase difference Δθ_(e)(ω_(m)) between voltage (V) and current (I) within the energy device from applying said mechanical excitations; and deducing the state-of-health of said energy device by comparing said electrically-induced phase difference Δθ_(e)(ω_(e)) with said mechanically-induced phase difference Δθ_(e)(ω_(m)).
 2. The method of claim 1, further comprising: taking a difference Δθ=Δθ_(e)(ω_(e))−Δθ_(e)(ω_(m)) between said electrically-induced phase difference Δθ_(e) and said mechanically-induced phase difference Δθ_(e)(ω_(m)); and using said difference as a basis to deduce said SoH.
 3. The method of claim 2, further comprising: deducing the state-of-health of said energy device by determining a magnitude of Δθ=Δθ_(e)(ω_(e))−Δθ_(e)(ω_(m)) in relation to 0°, wherein better state-of-health is deduced from said magnitudes being larger and poorer state-of-health is deduced from said magnitudes being smaller.
 4. The method of claim 2, further comprising: deducing that said energy device is approaching its end of life by determining that Δθ=Δθ_(e)(ω_(e))−Δθ_(e)(ω_(m))→0°.
 5. The method of claim 1, further comprising: at a baseline time, determining baseline-time phase differences Δθ_(e)(ω_(e)) and Δθ_(e)(ω_(m)) between current and voltage of said energy device; at a later time, determining later-time phase differences Δθ_(e)(ω_(e)) and Δθ_(e)(ω_(m)) between current and voltage of said energy device; and deducing said state of health by comparing said later-time phase difference with said baseline-time phase difference.
 6. The method of claim 5, further comprising: deducing said state of health by determining how much said later-time phase difference has diminished from said baseline time phase difference.
 7. The method of claim 6, further comprising deducing a magnitude of decline in said state of health substantially in proportion to said diminution.
 8. The method of claim 1, wherein: a difference between measurements of internal impedance Z_(INT) the energy device when said electrical excitation frequency ω_(e) is applied versus when said mechanical excitation frequency ω_(m) is applied, is determined according to: ${{{Z_{INT}^{e}\left( \omega_{e} \right)} - {Z_{INT}^{e}\left( \omega_{m} \right)}} = {\left( {\frac{\frac{1}{R_{el}\left( \omega_{e} \right)}}{\frac{1}{R_{el}^{2}\left( \omega_{e} \right)} + \left( {\omega_{e}{C_{b}\left( \omega_{e} \right)}} \right)^{2}} - \frac{\frac{1}{R_{el}\left( \omega_{m} \right)}}{\frac{1}{R_{el}^{2}\left( \omega_{m} \right)} + \left( {\omega_{e}{C_{b}\left( \omega_{m} \right)}} \right)^{2}}} \right) - {j\left( {\frac{\omega_{e}{C_{b}\left( \omega_{e} \right)}}{\frac{1}{R_{el}^{2}\left( \omega_{e} \right)} + \left( {\omega_{e}{C_{b}\left( \omega_{e} \right)}} \right)^{2}} - \frac{\omega_{e}{C_{b}\left( \omega_{m} \right)}}{\frac{1}{R_{el}^{2}\left( \omega_{m} \right)} + \left( {\omega_{e}{C_{b}\left( \omega_{m} \right)}} \right)^{2}}} \right)}}},$ where R_(el) designates an electrochemical resistance of said energy device and C_(b) designates a storage capacity of said energy device; further comprising: deducing said electrically-induced phase difference Δθ_(e)(ω_(e)) between said voltage V and said current I according to ${{Z_{INT}^{e}\left( \omega_{e} \right)} = \frac{V\left. \sqrt{}{\theta_{e}\left( \omega_{e} \right)} \right.}{I}};$ and deducing said mechanically-induced phase difference Δθ_(e)(ω_(m)) between said voltage V and said current I according to ${Z_{INT}^{e}\left( \omega_{m} \right)} = {\frac{V\left. \sqrt{}{\theta_{e}\left( \omega_{m} \right)} \right.}{I}.}$
 9. A apparatus for testing the state-of-health of an electrochemical energy device, comprising: an electrical excitation source for applying electrical excitations to the energy device at a predetermined electrical excitation frequency ω_(e); a mechanical excitation source for applying mechanical excitations to the energy device at a predetermined mechanical excitation frequency ω_(m); a measurement device for measuring an electrically-induced phase difference Δθ_(e)(ω_(e)) between voltage (V) and current (I) within the energy device from applying said electrical excitations; and said measurement device for further measuring a mechanically-induced phase difference Δθ_(e)(ω_(m)) between voltage (V) and current (I) within the energy device from applying said mechanical excitations; wherein the state-of-health of said energy device is deduced by comparing said electrically-induced phase difference Δθ_(e)(ω_(e)) with said mechanically-induced phase difference Δθ_(e)(ω_(m)).
 10. The apparatus of claim 9, wherein: the state-of-health of said energy device is deduced by taking a difference Δθ=Δθ_(e)(ω_(e))−Δθ_(e)(ω_(m)) between said electrically-induced phase difference Δθ_(e) and said mechanically-induced phase difference Δθ_(e)(ω_(m)).
 11. The apparatus of claim 10, wherein: the state-of-health of said energy device is deduced by determining a magnitude of Δθ=Δθ_(e)(ω_(e))−Δθ_(e)(ω_(m)) in relation to 0°, and wherein better state-of-health is deduced from said magnitudes being larger and poorer state-of-health is deduced from said magnitudes being smaller.
 12. The apparatus of claim 10, wherein: said energy device is approaching its end of life is deduced by determining that Δθ=Δθ_(e)(ω_(e))−Δθ_(e)(ω_(m))→0°.
 13. The apparatus of claim 9, wherein: baseline-time phase differences Δθ_(e)(ω_(e)) and Δθ_(e)(ω_(m)) between current and voltage of said energy device are determined at a baseline time; later-time phase differences Δθ_(e)(ω_(e)) and Δθ_(e)(ω_(m)) between current and voltage of said energy device are determined at a later time; said baseline-time phase difference is compared with said later-time phase difference to deduce said state-of-health.
 14. The apparatus of claim 13, wherein: said state of health is deduced by determining how much said later-time phase difference has diminished from said baseline time phase difference.
 15. The apparatus of claim 14, wherein a magnitude of decline in said state of health is deduced substantially in proportion to said diminution.
 16. The apparatus of claim 9, wherein: a difference between measurements by said measurement device of internal impedance Z_(INT) the energy device when said electrical excitation frequency ω_(e) is applied versus when said mechanical excitation frequency ω_(m) is applied, is determined according to: ${{{Z_{INT}^{e}\left( \omega_{e} \right)} - {Z_{INT}^{e}\left( \omega_{m} \right)}} = {\left( {\frac{\frac{1}{R_{el}\left( \omega_{e} \right)}}{\frac{1}{R_{el}^{2}\left( \omega_{e} \right)} + \left( {\omega_{e}{C_{b}\left( \omega_{e} \right)}} \right)^{2}} - \frac{\frac{1}{R_{el}\left( \omega_{m} \right)}}{\frac{1}{R_{el}^{2}\left( \omega_{m} \right)} + \left( {\omega_{e}{C_{b}\left( \omega_{m} \right)}} \right)^{2}}} \right) - {j\left( {\frac{\omega_{e}{C_{b}\left( \omega_{e} \right)}}{\frac{1}{R_{el}^{2}\left( \omega_{e} \right)} + \left( {\omega_{e}{C_{b}\left( \omega_{e} \right)}} \right)^{2}} - \frac{\omega_{e}{C_{b}\left( \omega_{m} \right)}}{\frac{1}{R_{el}^{2}\left( \omega_{m} \right)} + \left( {\omega_{e}{C_{b}\left( \omega_{m} \right)}} \right)^{2}}} \right)}}},$ where R_(el) designates an electrochemical resistance of said energy device and C_(b) designates a storage capacity of said energy device; said electrically-induced phase difference Δθ_(e)(ω_(e)) between said voltage V and said current I is deduced according to ${{Z_{INT}^{e}\left( \omega_{e} \right)} = \frac{V\left. \sqrt{}{\theta_{e}\left( \omega_{e} \right)} \right.}{I}};$ and said mechanically-induced phase difference Δθ_(e)(ω_(m)) between said voltage V and said current I is deduced according to ${Z_{INT}^{e}\left( \omega_{m} \right)} = {\frac{V\left. \sqrt{}{\theta_{e}\left( \omega_{m} \right)} \right.}{I}.}$ 